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SUMMARY: We present the radio detection of 18 BL Lac objects from our 
survey of over 575 deg of sky. These 18 objects are located within 20" of the 
X-ray position, of which 11 have a measured red-shift. All candidates are ra- 
dio emitters above — 1 mJy and fall within the range of existing samples on the 
two colour, OL ro vs CX ox , diagram with a transitional population of three (3) evi- 
dent. Two unusual sources have been identified, a candidate radio quiet BL Lac, 
RX J0140. 9-4130, and an extreme HBL, RX J0109. 9-4020, with Log (v pe ak) « 19.2. 
The BL Lac Log(N)— Log(S) relation is consistent with other samples and indicates 
the ROSAT All Sky Survey (RASS) could contain (2000±400) BL Lac objects. 
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1. INTRODUCTION 



BL Lac objects are a rare class of active galac- 
tic nuclei (AGN, see reviews by Kollgaard (1994), 
comprising only ~1.4% of the total AGN popula- 
tion (BL Lac known -1122, all AGN ~95971, Veron- 
Cetty & Veron 2006). Standard AGN unification 
models explain a BL Lac as a near line of sight, 
non-thermal (i.e. broadband synchrotron emission 
combined with inverse Compton scattering at higher 
frequencies), relativistic jet (i.e. Doppler boosting) 
produced by a FRI radio galaxy nucleus. Conse- 
quently, the optical spectrum dominated by contin- 
uum emission from the jet is characterised by weak 
emission or absorption lines. This makes the detec- 
tion of BL Lacs difficult at optical wavelengths using 



conventional AGN selection techniques. Most of the 
existing BL Lac samples have been initially selected 
at radio or X-ray wavelengths using large area sky 
surveys, followed by time consuming optical spec- 
troscopy to identify the BL Lac's in the sample. This 
led to the conclusion that BL Lac characteristics are 
bimodal, divided into radio selected BL Lacs (RSBL; 
i.e. the 1 Jy sample of Stickel et al. 1991) and X-ray 
selected BL Lac's (XSBL, i.e. the Einstein Extended 
Medium Sensitivity Survey (EMSS) sample of Rec- 
tor et al. 2001). 

A model proposed by Giommi & Padovani 
(1994) unifies BL Lac's as a single population, char- 
acterised by a spectral energy distribution (SED) 
with peak or cutoff of the synchrotron emission in the 
IR/optical (RSBLs) or in the UV/X-ray (XSBLs). 
To avoid confusion BL Lac studies use the classifi- 
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cation suggested by Padovani & Giommi (1995) of 
low energy peak BL Lacs (LBLs, mostly RSBL) or 
high energy peaked BL Lacs (HBLs, mostly XSBL). 
Both classes exhibit similar extreme "Blazar" activ- 
ity (see Urry & Padovani 1995) — rapid, irregular 
variability, high optical polarization, radio core dom- 
inated morphology (Laurent-Muehleisen et al. 1993), 
flat radio spectra, super-luminal motion and emis- 
sion from the radio to gamma ray bands (Sam- 
bruna et al. 1996). Given the EMSS and 1 Jy sam- 
ples represent extremes of the BL Lac population 
Giommi & Padovani (1994) predicted the existence 
of BL Lacs with properties intermediate to those of 
LBLs and HBLs — termed intermediate BL Lacs 
(IBLs). 

Currently our understanding of BL Lac ob- 
jects is far from complete, due mainly to their 
rarity limiting the size of statistically complete 
samples. One area that is poorly understood is 
the luminosity evolutionary behavior of BL Lacs. 
While other classes of AGN show a strong trend 
of positive evolution, BL Lacs do not. LBLs 
and IBLs exhibit no or a slight positive evolution 
(Stickel et al. 1991; Rector & Stocke 2001; Laurent- 
Muehleisen et al. 1999) while HBLs exhibit a neg- 
ative evolution (Rector et al. 2001), meaning they 
are less luminous/frequent in the past than now. Al- 
though the properties of HBLs and LBLs are dif- 
ferent, they are not sufficient to determine if the 
properties reflect two different populations of AGN's 
or the extremes of a single continuous AGN pop- 
ulation. The role of IBLs remain uncertain (Rec- 
tor et al. 2003) as their observational properties are 
not well understood due to the small number known 
(Laurent-Muehleisen et al. 1999). 

In the last 18 years (since the launch of 
ROSAT), the number of BL Lacs has doubled ev- 
ery six to seven years (Fig. 1) ("IS M direct result of 
efficient selection techniques (Schachter et al. 1993; 
Wolter et al. 1997). Stocke et al. (1989) noted that 
since BL Lacs are radio and X-ray loud, their broad- 
band spectral properties could be used to select out 
candidate BL Lac objects. In a purely X-ray survey, 
efficiencies of <5% are reached (Stocke et al. 1991), 
increasing to ^10% at Fx>10~ 12 ergs cm~ 2 s _1 
(Perlman et al. 1996). In a purely radio survey ef- 
ficiencies of ~7% are reached (Stickel et al. 1991), 
increasing to about 12% if flat spectrum sources 
are pre-selected. In a combined radio and X-ray 
survey efficiencies of 10% to 20% or higher can be 
reached (Schachter et al. 1993; Wolter et al. 1997). 
Recently, positional correlation of radio and X-ray 
catalogues has been undertaken: Deep X-ray Ra- 
dio Blazar Survey (Perlman et al. 1998), the REX 
Survey (Maccacaro et al. 1998) and the ROSAT All 
Sky Survey - Green bank Survey (RGB; Laurent- 
Muehleisen et al. 1999). A large sample of 501 radio- 
selected BL Lac candidates has been compiled by 
Plotkin et al. (2008) using Faint Images of the Ra- 
dio Sky at Twenty- Centimetres (FIRST) radio sur- 
vey in combination with the Sloan Digital Sky Survey 



(SDSS) optical spectroscopy. 
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Fig. 1. Plot of the number of published BL Lac in 
the literature as a function of time, based on Veron- 
Cetty & Veron (1981 1985, 1987, 1989, 1991, 1993, 
1996, 1998, 2000, 2001, 2003, 2006). For 2006 
— 1.2% of all known AGNs are BL Lac, compared 
to -88. 8% for QSO and -10.0% for Syl. As new 
ROSAT BL Lacs are discovery the number of known 
BL Lacs should exceed -2000 by 2010. 



Potentially, the largest complete catalogue 
of BL Lac objects remains mostly undiscovered in 
the ROSAT All Sky Survey (RASS), a fact first 
noted by Stocke et al. (1989). Various studies 
have begun the task of identifying RASS samples of 
BL Lacs, the RGB sample (127 BL Lacs) of Laurent- 
Muehleisen et al. (1999) and SDSS/RASS sample of 
Anderson et al. (2007) with 200 new BL Lacs. The 
limiting factor for RASS sample sizes is the time con- 
suming task of optical identification with a 15-20" 
position error. 

Presented in this paper is a sample of 18 new 
RASS BL Lacs (including five (5) candidate objects 
and five (5) objects verified from the literature) dis- 
covered as part of a radio detection survey and ESO 
optical identification program for a complete sub- 
sample of 695 RASS sources. Section 2 describes the 
RASS sub-sample selection, radio follow-up program 
and details how optical identifications were made and 
the criteria used to select objects as BL Lacs. In 
Section 3 we describe the sample, including a table 
of BL Lac properties. In Section 4 we examine the 
properties of the sample and compare them to sam- 
ples from the literature. Finally, in Section 5 we 
discuss the findings with Section 6 forming the con- 
clusion. 
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2. THE RASS SUB-SAMPLE 



The ROSAT X-ray telescope completed the 
first All Sky Survey with an imaging X-ray tele- 
scope in 1991, resulting in two RASS Catalogues 
(unpublished RASS revision one catalogue, ^54000 
sources, see Voges (1993) and various RASS revision 
two catalogues, ^160 000 sources, see Voges (1997) 
and Voges et al. (1999, 2000)). To characterise the 
source population being detected a complete sub- 
sample of four (4) Key Fields, Table 1, was selected 
(|b| >-25°, <5<10°, 695 sources, covering 575 square 
degrees) from the revision one catalogue for optical 
identification (Danziger et al. 1990; European South- 
ern Observatory Key Field Program). To improve 
the multi-wavelength coverage and to assist the op- 
tical program all Key Field sources where observed 
at 4.8 GHz (~1 mJy limit and <~1" positional accu- 
racy). 

Table 1. Summary of the RASS Key Field proper- 
ties. 



Key Field 


RA 


Dec 


Number of Sources 


I 


01 h 


-40° 


261 


II 


05 h 


-57° 


227 


III 


n h 


-27° 


71 


IV 


13 h 


+00° 


134 



The optical and radio follow-up identified a 
complete sample of 18 BL Lacs using X-ray posi- 
tions from the unpublished RASS revision one (1) 
catalogue. A positional correlation with the RASS 
revision two (2) catalogue establish all of the objects 
are real detections. 



2.1 Radio follow-up 



Radio observations of the BL Lac sample 
(695) were undertaken as part of the 4.8 GHz 
radio detection program (Anderson 2003; Ander- 
son et al. 1994). Briefly, the Australia Telescope 
Compact Array (ATCA) observed the RASS Key 
Field sources in fields I, II and III, while the Very 
Large Array (VLA) observed sources in field IV. 
For the ATCA all observations used a 1.5 km array 
configurations plus antenna 6 and were made using 
a 128 MHz bandwidth divided into 32 channels of 
width 4 MHz. Closest match for the VLA is the C 
array configuration, with a bandwidth of 50 MHz. 
All radio detections are located within 2' of the 
pointing center where bandwidth smearing effects 
are extremely small. The primary flux calibrators, 
PKS 1934-638 (ATCA) and 3C286 (VLA), are tied 
to the VLA absolute flux scale of Baars et al. (1977). 

The radio data has been reduced using stan- 



dard reduction procedures within the Astronomical 
Image Processing Software (AIPS) package. The ob- 
servations reached a five sigma detection limit of 
-2 mJy. All 4.8 GHz flux densities listed in Table 2 
have been primary beam corrected. 

The radio positions were correlated with the 
NRAO VLA Sky Survey (NVSS) radio survey (Con- 
don et al. 1998) to determine 1.4 GHz flux densi- 
ties and errors for BL Lacs north of a declination of 
-42°. Ten coincidences were found, with flux den- 
sities listed in Table 2. The efficiency of combined 
radio (NVSS) and X-ray (RASS) can be as large as 
35% (Beckmann et al. 2003). A similar compari- 
son with the Sydney University Molonglo Sky Survey 
(SUMSS; Mauch et al. 2003) to determine 843 MHz 
flux densities found four source, listed in Table 2. 

Four BL Lacs were verified by a positional cor- 
relation of the Key Field radio positions with BL Lac 
objects in the NASA/IPAC Extragalactic Database 
(NED). 



2.2 Optical follow-up and identification 



Preliminary optical candidates were selected 
using optical finder charts extracted from the COS- 
MOS catalogue (Yentis et al. 1992). Candidates were 
prioritised according to distance from the RASS po- 
sition and optical magnitude, with priority given to 
optical sources associated with a radio counterpart. 
Optical spectroscopy was then used to identify one 
or more optical candidates until the most likely iden- 
tification was established. 

X-ray selected BL Lacs are X-ray loud and ra- 
dio loud (Stocke et al. 1991). In practical terms this 
means all EMSS (Rector et al. 2001) and ESS (Perl- 
man et al. 1996) BL Lacs have a 4.8 GHz flux den- 
sity brighter than 1 mJy. Hence all RASS BL Lacs 
should be positionally associated with a radio coun- 
terpart brighter than 1 mJy. If all BL Lacs are 
required to be associated with a radio source then 
false identification from stellar source with feature- 
less spectra would also be minimised. In the Key 
Field sample only one stellar X-ray sources (Ander- 
son& Filipovic 2007) was detected as a radio emitter. 

BL Lac object classification was based on 
the technique developed for the Einstein Extended 
Medium Sensitivity Survey (Stocke et al. 1991). A 
source was classified as a BL Lac if it had: 

(i) A featureless optical spectrum i.e. no observed 
emission lines of equivalent width >0.5 nm. 

(ii) A Ca II break of less than 25% to distinguish 
between a non-thermal AGN spectrum and 
the stellar spectrum of an elliptical galaxy. 

(iii) A 4.8 GHz flux density greater than 1 mJy. 

Any candidate BL-Lac object was required to 
meet criteria three and have at-least a featureless 
spectrum. 
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More details about the optical program is de- 
scribed by Danzigcr ct al. (1990). We recognise that 
the strength of the Calcium break is depending on 
the luminosity of the object: the strength of the 
Calcium break is anti-correlated with the luminosity 
(Landt et al. 2002; Beckmann et al. 2003). Hence, 
some BL Lac objects could exhibit Calcium breaks 
exceeding 25% (Anderson et al. 2007). In addition, 
blazars can show emission lines. As the BL Lac ob- 
jects are a subclass of the blazar population and as 
there is apparently a smooth transition from BL Lacs 
to Flat Spectrum Radio Quasars (FSRQ), one has to 
be very careful when applying strict constraints on 
the optical spectra. While we acknowledge that we 
might have missed border-line BL Lac objects, the 
use of criteria three should help reduce the number 
missed. 



3. THE BL LAC SAMPLE 



A total of 18 BL Lacs have been identified, 
their characteristics are summarised in Table 2; 
Column (1) is the table index for cross referenc- 
ing with the index in Table 4; Column (2) is the 
ROSAT All Sky Survey Name; Column (3) is the 
COSMOS J band magnitude; Column (4) is the red- 
shift; Column (5) is the radio position right ascen- 
sion for J2000; Column (6) is the radio position dec- 
lination for J2000; Column (7) is the ATCA 4.8 GHz 
flux density and error in mJy; Column (8) is the 
1.4 GHz flux density in mJy from either the Parkes 
90 catalogue or the NVSS (only available for sources 
north of a declination of -42°) and associated errors; 
Column (9) is the reference flag defined at the end 
of Table 2 and indicates if a source is a candidate or 
confirmed Key Field BL Lac or is from the literature; 
Column (10) is the RASS count rate and associated 
error and column (11) is the Total Hydrogen Column 
Density in units of 10 20 atoms cm -2 . For additional 
details consult the source notes at the end of this 
section. 

A total hydrogen column density, Njj, value 
(Col. 11) is needed to correct the calculated ROSAT 
Position Sensitive Proportional Counter (PSPC) flux 
for galactic absorption. Values were estimated for 
Key Field sources north of a declination of -42° us- 
ing the Bell Labs Hi survey of Stark et al. (1992). 
For Key Field sources south of a declination of 
42°, estimates were taken from the Parkes Hi sur- 
vey of Cleary et al. (1979). Only values from 
Stark ct al. (1992) are not affected by stray radia- 
tion entering the far field antenna pattern. 

Thirteen (13) BL Lacs, including five (5) 
candidates, have been optically identified by 
Danziger et al. (1990) and are listed in Table 2. Up to 
June 2009, a positional correlation of the radio posi- 
tions with objects in the NASA/IPAC Extragalactic 
Database (NED) identified five (5) additional sources 



classified in the literature as BL Lacs and listed in 
Table 2. 

Multi-epoch 4.8 GHz radio observations, 
based on a set of 40 sources, indicates the flux den- 
sity error is well represented by or = O.OSS^ 1 +0.48, 
with Sr the flux density in mJy over the range of 
1.5 mJy to 1 Jy. Above 30 mJy the flux error is con- 
stant at about 8%, rising to 40% at 1.5 mJy, with 
errors for individual flux densities given in Table 2. 

The number of randomly associated radio 
counterparts for the entire Key Field sample of 695 
RASS objects is estimated at less than two, based on 
4.8 GHz radio source counts. This makes the chances 
of a BL Lac being a random association extremely 
small. 

The Key Field sample is complete to approx- 
imately (9 ± 2) x 10 -13 ergcm~ 2 s _1 , at which, the 
prior counts of Wolter et al. (1991) predict a sky 
density of 0.030±0.008 per degree squared or 17±5 
BL Lacs over the 575 deg 2 of sky covered by the 
four Key Fields. This value compares well to the ob- 
served value of 18±4 BL Lacs observed, indicating a 
high level of completeness. The number of BL Lacs 
missed is estimated to be less than two. 



3.1. Positional accuracy 



Astrometric accuracy of the radio and X-ray 
positions determined by Anderson (2003) are 1.1" 
and 30" respectively. For the radio positions the 
accuracy was determine by comparing the positions 
of a standard set of 40 sources observed at two or 
more epochs, giving a deviation of A6 — 0.8" and 
Aacos(5) = 0.7". RASS positional accuracy was de- 
termined by comparing the X-ray and radio positions 
for the ~ 120 radio emitting X-ray sources identified 
by the Key Field program, with 95% of the radio 
sources located within 40" of the X-ray position. 

Comparison of the difference between the 
RASS and ATCA radio positions for the 18 BL Lacs 
in Fig. 2 indicates 90% of identified sources are lo- 
cated within 15". This value is smaller than the 30" 
error for the entire Key Field sample and is probably 
due to BL Lac morphology being compact. 
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Fig. 2. The positional difference between the RASS 
X-ray and the ATCA radio position. All BL Lac's 
are located within 20" from the X-ray position. 



to be a BL Lac object. Veron-Cetty & Veron 
(2006) confirmed RX J2324.7-4041 (also known as 
1ES 2322-409) as a BL Lac. 

Xue et al. (2000) investigated RX J2319. 1-4206, 
showing it to be a BL Lac with a complex radio mor- 
phology and it is also known as PKS B2316-423 or 
PMN J2319-4206. The compact core has a ATCA 
4.8 GHz flux density of 207 mJy, while the extended 
lobe has a 237 mJy flux density, giving a total of 
444 mJy that compares well to the PMN single dish 
value of 595 mJy. We also note that this object is 
very close to the galaxy cluster Abell SI 111 located 
1.1' away. This would be indeed very interesting, 
as BL Lac objects usually seem to avoid cluster 
environments (Wurtz et al. 1993, 1997). However, 
upon closers inspection we found that the redshift of 
RX J2319.1-4206isz=0.055and that of Abell Sllll 
is z=0.045, which gives a difference in distance of 
~50 Mpc. Thus, this BL Lac is not likely inside the 
galaxy cluster. 



4. PROPERTIES OF OUR SAMPLE 



3.2. Source notes 



BL Lacs classified as candidate objects 
include the RASS source RX J0049.7-4151, 
RX J0045.8-3158, RX J0059. 5-3510, RX J0054.8-2455 
and RX J1357.2-0146. 

Objects confirmed as BL Lacs by 
the ESO optical identification program in- 
clude RX J0109.9-4020, RX J0140.9-4130, 
RX J0543.9-5532, RX J0506.9-5435 (also known as 
1ES 0505-546), RX J0528.8-5920 (located next to a 
field galaxy), RX J0528.8-5920, RX J0439.0-5558, 
RX J1357.6+0128 and RX J1103.6-2329. Multi 
epoch radio data for RX J0439. 0-5558 reveals a 
variability of 28% or 3 mJy. RX J0109.9-4020 
was independently confirmed as a BL Lac 
(see Schwope et al. 2000), with the radio 
source PMN J0110-4020 one arcminute away. 
RX J0543. 9-5532 was independently identified in 
NED as a BL Lac (see Fischer et al. 1998). 
RX Jl 103.6-2329 was independently confirmed 
by Laurent-Muehleisen et al. (1993) and Perl- 
man et al. (1996) as a BL Lac, and is also known 
as 1ES 1101-232, and PMN J1101-2329 which has 
a 4.8 GHz flux density of 66±11 mJy that compares 
well to the ATCA value. 

Five additional objects are classified as 
BL Lacs in the literature. Bauer et al. (2000) identi- 
fied RX J0040. 3-2719 as a BL Lac and it is also coin- 
cident with the radio source NVSS J004016-271912. 
RX J0043.4-2639 was confirmed as BL Lac by Cris- 
tiani et al. (1995), also known as PMN J0043-2639 
which has a slightly higher 4.8 GHz flux density of 
81±11 mJy compared to the ATCA value 68±6 mJy. 
Bade et al. (1994) confirmed RX J1057.8-2753 



The characteristics of the sample are ex- 
amined by comparing it with the Einstein EMSS 
BL Lac sample (Gioia et al. 1990; Rector et al. 2001), 
the 1 Jy BL Lac sample (Stickel et al. 1991; Rec- 
tor & Stocke 2001) and the Einstein Slew survey 
(ESS) BL Lac sample (Perlman et al. 1996). The 
1 Jy sample is taken to be representative of the prop- 
erties of Low Energy Peaked BL Lacs (LBLs) and 
the EMSS of the properties of high energy peaked 
BL Lacs (HBLs). The ESS exhibits properties that 
are intermediate between LBLs and HBLs. Compar- 
ison was made of the distributions of v pea k (peak fre- 
quency of the synchrotron emission), redshift, a no 
versus aox, median spectral index and the X-ray 
BL Lac Log(N)-Log(S). The median values for the 
four samples are given in Table 3. 



Table 3. Median properties of the four BL Lac sam- 
ples. The median redshift value for this survey is an 
estimate as only 55% of sources have redshifts. 





EMSS 


1 Jy 
Radio 


ESS 


This 
Study 


Redshift 


0.30 


0.55 


0.16 


0.31 


Log (Sx/Sr) 


-4.8 


-6.9 


-4.6 


-5.0 


Log(u peak ) 


16.2 


13.8 


15.4 


15.9 



4.1. Peak frequency of the synchrotron emis- 
sion 



The turnover frequency at which the spec- 
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tral energy distribution (SED) peaks can be used 
to distinguish LBLs from HBLs. Following Sam- 
bruna et al. (1996) and Landau et al. (1986) the key 
field sample flux densities for the radio (Sjj) optical 
(So) and X-ray (Sx) were fitted with a logarithmic 
parabola of the form: 

Log {yS v ) = a(Log (v)) 2 + bLog [y) + c (1) 

with v the observed frequency. The peak frequency 
is then given by Log (is pea k) = —b/2a and was k- 
corrected to the rest frame using v rest = v (1 + z) 
were v rest is the rest frame frequency. 

For the seven (7) BL Lacs with no measured 
redshift, a value equal to the sample median (z=0.31, 
Table 3) is used. While a median will likely under- 
estimate the true redshift (i.e. nearby BL lacs (z< 
0.2) tend to have measured redshift due to the host 
galaxy being resolved, while objects with no measur- 
able redshift tend to be at z<l) any associated error 
in ctRx due to a K-correction would be smaller then 
0.04 for z<l. 

Only radio, optical and X-ray data from Ta- 
ble 2 was used. Flux densities were calculated for the 
PSPC count rate at 2 keV and for the COSMOS op- 
tical magnitude at 641 THz using S =3.95x 10~ 23 
lcr o.4m(j) adopted from Fukugita et al. (1995) for 
So values calculated within the NED database. An 
example of the SED and functional fit for one source 
is given in Fig. 3. The k-corrected peak values are 
listed in Table 4. A fit to source 7 (RX J0109.9-4020) 
was problematic, giving a very high peak frequency 
that has been placed in brackets in Table 4. Ex- 
amination of the broad band fluxes shows the X-ray 
count rate is abnormally high (see Sect. 5). 
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Fig. 3. Example of the logarithmic parabola fit- 
ted ( equation 1) to the spectral energy distribution of 
source RX ,100^3.^-2639 , Table 4, with the observed 
Log(v peak ) =14-81. 



Fig. 4 compares the two point radio to X-ray 
spectral index, a R x, versus Log(i/ pea fe) diagram for 



this sample (Table 4) with that of the one (1) Jy 
Radio and the EMSS using values from the liter- 
ature (Sambruna et al. 1996). Clearly aux and 
\jOg{v pea k) are correlated. A linear regression analy- 
sis of the Key Field data, excluding one anomalous 
point (source 7 marked with \ in Table 2 and plotted 
as solid diamonds in Fig. 4), results in: 

ctRx = -0.091 x Log {v pea k) + 2.12 

In good agreement to the fit of Dong et al. (2002), 

a RX = -0.09535 x Log (v peak ) + 2.16065, 

based on the EMSS and 1 Jy radio samples. 

Two lines defined by aux=0-75 and 
Log(i'p eQ fe)=14.7 in Fig. 4 can be used to delineate 
four regions. HBLs occupy the lower right quadrant 
and LBLs the upper left quadrant. This distinc- 
tion was first suggested by Dong et al. (2002) and 
is equivalent to using Log(Sr/Sx)=-5.5 (Laurent - 
Muehleisen et al. 1999) as a criteria. Intermediate 
BL Lacs occupy the region defined by a R x—0-7 to 
a RX =0-8. 

The Key Field sample is composed mostly of 
13 HBLs, three IBLs and two others falling within 
the region for LBLs. A median Log(i/ pe afc)=15.9 
(Table 3) is close to the EMSS value of 16.2, both 
are typical of HBLs. As expected the ESS with 
IBLs has a median value \jOg(v pea k) = \^ .4 that falls 
close to halfway between the LBLs (1 Jy sam- 
ple) Log(i^ eQfe ) = 13.8 and the HBLs Log(i/ peafc )=16.2 
(EMSS). 



Table 4. Values of the spectral indices for radio 
to optical, optical to X— ray and radio to X— ray, the 
logarithm of the frequency for the SED turnover peak 
and the type of BL Lac, LBL, IBL and HBL. Note 
that Col. 1 (Index) is for cross referencing with Ta- 
bic 1. 



Index 


OlRO 


CtOX 


OLRX 


Log(^ peafc ) 


Type 


1 


0.42 


1.12 


0.70 


15.9 


HBL 


2 


0.41 


1.52 


0.83 


14.8 


LBL 


3 


0.38 


0.99 


0.63 


16.6 


HBL 


4 


0.42 


1.09 


0.69 


16.0 


HBL 


5 


0.35 


1.22 


0.68 


15.5 


HBL 


6 


0.50 


1.03 


0.72 


16.4 


IBL 


7 


0.56 


0.66 


0.64 


(19.2) 


HBL 


8 
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4.2. Redshift distribution 

The Key Field sample has 11 BL Lacs with 
measured redshifts (eight (8) have z<0.5 and three 
(3) have z>0.5) and seven with no redshift. Six of 
the seven with no redshift arc bright (i.c mj<18) 
and have Log(S m (j) / S x) typical of objects in the 
EMSS with a z<0.8 (Stocke et al. 1991). The me- 
dian redshift of the sample (z=0.31) is consistent 
with the EMSS (z=0.30) and is placed about halfway 
between the median value for the ESS (z=0.16) and 
1 Jy (z=0.55) samples. Similarly the median z and 
Log(Sjf /Sfi) for the three samples and the Key Field 



value of Log(Sjf /Sfl)=-5.0 corresponds to an inter- 
polated median z of 0.35±0.4. 

The Key Field sample indicates a sample of 
HBLs has a systematically lower redshifts than a 
LBLs (i.e. 1 Jy) sample. This finding is consistent 
with the unification model of Fossati et al. (1997), 
which predicts HBLs to have systematically lower 
redshifts compared to LBLs. Meaning that HBLs 
have lower intrinsic luminosities than LBLs. A 
slight correlation was found for the Key Field and 
EMSS sample between z and Log(^ pea fe), indicating 
~Log{v pea k) is increasing with z. The same variables 
do not correlate for the 1 Jy sample. 




Fig. 4. The positional difference between the RASS X-ray position and the ATCA radio position. All 
BL Lac 's are located within 20" from the X-ray position. 
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Fig. 5. Radio to optical versus optical to X-ray spectral index plot for BL Lac objects. The Key Field 
sample, including candidate objects, are plotted as a large diamond and triangle, respectively. Other samples 
shown are for comparison and are drawn from X-ray selected samples; Einstein EMSS (Giommi et al. 1995), 
Einstein Slew survey (ESS; Perlman et al. 1996), other ROSAT BL Lacs and radio selected BL Lacs; 1 Jy 
and S5 (Stickel et al. 1991). The regions defining different types of objects, and the definition of Uro and 
aox, are taken from the Einstein EMSS (Stocke et al. 1991). A one magnitude optical variation is indicated 
in the top right corner. 



4.3. Spectral index 



The method of Brinkmannfe Siebert (1994) 
was used to estimate a spectral index for each 
BL Lac using hardness ratios. The hardness ratios 
for each BL Lac have large errors that make the 
estimated spectral index values uncertain. There- 
fore, only the sample median value will be pre- 
sented. For the 18 sources the median energy index 
is 1.42±0.80 and compares well to values from the 
literature, 1.32±0.40 for a RASS sample (Laurent- 
Muehleisen et al. 1999) and 1.47±0.40 for the EMSS 
(Rector et al. 2001). 



4.4. otRo versus aox diagram 



Fig. 5 is the uro versus aox plot 
for the Key Field and comparison sam- 
ples. The two point spectral index is calcu- 
lated using a_ RO =Log(Si ? /S )/Log(i/_ R /^ ) and 



aox=Log(So/Sx)/Log(i'o/ l/ x) with flux densities 
calculated at v R =A.S GHz, ^ =641 THz (COSMOS 
Ill-aJ Photographic B magnitude is ~ B A=4400 A, 
see Sect. 4.1 for further details) and vx = 5.40xl0 17 
Hz (X-ray flux density is calculated for 2 keV). The 
flux densities have been k corrected to the rest frame 
(i.e. S rest =Sx (l+z) Q_1 ), assuming a/{=0.0, ao = l-0 
and ax=l-2), with the median sample redshift (Ta- 
ble 3) used if none are available. The final values 
are listed in Table 4, along with a^x- Comparison 
values are based on data in the literature. 

The one anomalous Vp ea k 

value from Sect. 4.1 
corresponds to an extreme outlier in the apto versus 
aox plane above the main BL Lac region. The value 
is flagged in Table 4 by brackets and is likely unrep- 
resentative of the overall distribution (see Sect. 5). 

The other 15-points (4 candidate and 11 con- 
firmed), likely HBLs, are positioned within the re- 
gions occupied by the EMSS and the 1 Jy samples. 
Two points, likely IBLs, are located in the radio quiet 
region, but are still within the region occupied by 
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ESS BL Lacs. 



4.5. X-ray Log (N) - Log (S) 



The X-ray source counts (Fig. 6) are consis- 
tent with the EMSS and the ESS counts. They have 
a slight bump at fx = 10~ n erg cm~ 2 s _1 reported 
by Maccacaro et al. (1989) and Nass et al. (1996). 
The EMSS counts fall well below the Key Field 
level as does the ESS count. This is partially 
due to incompleteness of the EMSS to bright X-ray 
sources and the ESS to faint X-ray sources between 
^2xl0~ 12 erg cm -2 s -1 and 10 -11 erg cm -2 s _1 . 



0.1 



5? 

"U 



X 
LL 



0.01 



£. 0.001 




x EMSS 
+ ESS 
o This study 



0.0001 



10-13 10-12 10-11 10-1° 
F x (erg cm -2 s _1 ) 



10-09 



Fig. 6. X-ray Log (N)-Log(S) for the Key Field 
BL Lacs and comparison samples. All error bars are 
based on the count error. 



Taking the flux limit of the RASS as 
~3xl0™ 13 erg cm~ 2 s _1 (0.1-2.4 keV) corresponds 
to a sky density of 0.04 BL Lacs per square de- 
gree. At this density the RASS will contain 
(2.6±0.4)xl0 3 BL Lacs for the whole sky, reducing 
to (2.0±0.4)xl0 3 at high galactic latitudes. 



DISCUSSION 



No indication of variability was found, while using 
only fluxes close to the X-ray flux epoch (1991) does 
not alter the Log(y pea k) value. The only notice- 
able difference is the RASS flux is significantly high 
relative to the other optical and X-ray fluxes i.e. 
c*ox = 0.66 while aox > 1 for the 17 other ob- 
jects. Alternatively, RX J0109. 9-4020 could be an 
example of a new population of BL Lacs, possibly a 
ultra high energy peaked BL lac predicted by Ghis- 
ellini (1999) to have a peak frequency v pea k located 
at higher frequencies than HBLs, Log(y„ ea ^) > 18. 
Nieppola et al. (2006) and Wu et al. (2009) showed 
for a large BL Lac sample their is a population of 
extreme HBL with Log(v pea k) values extending out 
to 21.5. 

Traditionally the ano-dox plot has been 
used to demonstrate the divisions of the popula- 
tion. These boundaries have proven successful in 
identifying new BL Lac objects, achieving efficien- 
cies of ~20% (e.g. Einstein Slew Survey (ESS); Perl- 
man et al. 1996; Schachter et al. 1993). If the EMSS 
colour-colour boundaries had been strictly applied 
to the Key Field BL Lacs then up to seven BL Lacs 
would have been missed. 

A problem facing BL Lac samples is the use 
of selection criteria and the associated biases in- 
troduced. An example of this is the observed bi- 
modality between LBLs (1 Jy) and HBLs (Key Field 
and EMSS) which is partly due to BL Lac selec- 
tion, in that the 1 Jy and EMSS/Key-field samples 
contain objects biased to the selection of BL Lacs 
using radio or X-ray properties only. Laurcnt- 
Muehleisen et al. (1999) showed the ROSAT Green 
bank BL Lac sample has no bimodality owing to a 
large number of intermediate BL Lacs. 

By cutting off the Key Field sample to sim- 
ulate an X-ray and Radio selected sample should 
reproduced the observed bimodality. The EMSS 
limit in the ROSAT band is about 3 x 10~ 2 ^Jy 
(Laurent-Muehleisen et al. 1999), giving a median 
Log(Sx/S_R)=-4.8, which is exactly the same as the 
EMSS (HBL) value, Table 3. The only bright radio 
source in the sample was taken to be representative 
of the 1 Jy sample, giving Log(Sx/S_R)= -6.7 for the 
total flux, which is the same as the 1 Jy (LBL) me- 
dian value of -6.9. Clearly even the Key Field sample 
can reproduce the observed bimodality between HBL 
and LBL. 



A candidate radio quiet BL Lac, RX J0140.9- 
4130, with a 1.8 mJy flux density occupies the radio 
quiet AGN region of the olro-clox plot at ano=0A5 
(Fig. 5). A total of 16 BL lacs in the sample have 
ctB,o > 0.3, more than double the 0.15 value. This 
contrasts to most X-ray selected samples that tend 
to be X-ray loud and radio loud (Stocks et al. 1991). 

The unusual source, RX J0109.9-4020, has a 
SED with a very high synchrotron peak frequency 
with Log (i/peak) ~ 19.2. Using all available flux 
measurements we looked for possible abnormal flux 
values that might explain the high Log (v pe ak) value. 
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CONCLUSION 



The primary result of this paper is the sample 
of 18 BL Lacs presented in Table 2. Based on the 
sample properties (Table 4) their are two LBL, three 
IBLs and 13 HBLs, making HBLs over represented 
compared to LBLs that make up only 11% of objects, 
as expected for a primarily X-ray selected sample. 
By applying appropriate X-ray and radio flux limits 
to the sample we where able to reproduce the ob- 
served bimodality between LBLs (1 Jy) and HBLs 
(EMSS). Two unusual source have been identified, 
a candidate radio quiet BL Lac, RX J0140. 9-4130, 
and an extreme HBL, RX J0109.9-4020, with a very 
high synchrotron peak frequency Log (v P eak) ~ 19.2. 
A slight positive correlation was found between fpeak 
and redshift. 

A correlation between anx and 
Log^peafc) was also found defined by the fit 
anx = — 0.091Log (V pea fc) + 2.12 and is in good 
agreement to the fit of Dong et al. (2002). The 
fit can be used to directly calculate v pe ak when no 
optical flux is available. 

Two main problems are faced by BL Lac stud- 
ies, understanding the SED and how it relates to 
the properties of the underlying source population, 
and the small number of known objects. A number 
of major programs are constructing large BL Lac 
samples (w 500) by combining the RASS (Ander- 
son et al. 2007) or FIRST (Plotkin et al. 2008) sur- 
veys with the SDSS for identification. At a flux 
limit of ~3xl0~ 13 erg cm" 2 s" 1 (0.1-2.4 keV) the 
RASS contains ^2592 BL Lacs for the whole sky, 
reducing to ^2000 at high galactic latitudes. Even 
a complete sample of only 500 objects would al- 
low major progress to be made in understanding 
the BL Lac phenomena. Our concept of existing 
BL Lac objects based on radio or X-ray selection 
may need to be modified to fit into potentially new 
BL Lac populations suggested by optical selections 
Londish et al. (2002) and Collinge et al. (2005). 
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Table 2. New BL Lacs identified by Key Field program. Table references for column nine: (1) candidate Key Field BL Lac; 
(2) Confirmed Key Field BL Lac; (3) Confirmed BL Lac from the literature; (4) Schwope et al. (2000); (5) Fischer et al. (1998); 
(6) Laurent-Muehleisen et al. (1993) and Perlman et al. (1996); (7) Bauer et al. (2000); (8) Cristiani et al. (1995); (9) 
Xue et al. (2000); (10) Bade et al. (1994); (11) Veron-Cetty & Veron (2006). 
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f RX J2319. 1-4206 has a ATCA 4.8 GHz flux density of 207 mjy for the core and 237 mjy for the lob, 444 mjy in total. 
RX J0109.9-4020 1.4 GHz flux density is for ATESP J010956-402051 from Prandoni et al. (2000). RX J0439.0-5558 has a 
843 MHz flux density of 21.2+1.3 mjy for SUMSS J043901-555840 from Mauch et al. (2003), and includes RX J0506. 9-5435 
with 17.1+1.1 mjy, RX J0528. 8-5920 with 18.8+1.0 mjy, RX J0543. 9-5532 with 41+2 mjy. 
*For more details see Section 4.1. 



M.W.B. ANDERSON and M.D. FILIPOVIC 



PAJIMO HETEKIIHJA 18 RASS BL LAC OBJEKATA 



M. W. B. Anderson 1 2 and M. D. Filipovic 2 



1 Sydney Observatory, PO Box K346, Haymarket, Sydney, NSW 1238, Australia 

2 University of Western Sydney, Locked Bag 1797 
Penrith South, DC, NSW 1797, Australia 

yuK ... 

Opuzunajiuu uayunu pad 

y OBoj CTy^HjH npcucTaBJbaMO pa,HHO /je- 
TCKrpijy 18 hobhx BL Lac oSjeKaia ca Hainer 
nperjie^a He6a or npeKO 575 kb a,np sthmx ctc- 
ncHH. Obhx 18 oSjeKaTa cy jion,npaHH y okojihhh 
oa 20" oa X-ray no3niiHje. 3a 11 oa 18 06- 
jeKaTa hms-mo H3MepeHn npBeHH noiiaK. Cbhx 18 
KaH/rH^aTa cy pa^no ejviHTcpn H3K&A ~1 mJy h 
noAnaAajy y no3HaTy rpynairnjy oa u Ave 6oje", 
a ro vs a ox , AnarpaM ca gbh^ghthom TpaH3n- 



ithohom nonyjiaiiHjoM 3. JleTeKTOBanH cmo h 
ABa BeoMa HeoSnMHa n3Bopa: KaH/rj/maT 3a pa- 
Avto THxy BL Lac, RX J0140. 9-4130, h ckctpgmhh 
XBJI, RX J0109. 9-4020, ca Log{v peak ) w 19.2. 
BL Lac Log(N)-Log(S) pejiairnja je KOHCHCTeHTHa 
ca ocTajiHM y3oprj;HMa h HH^HiiHpa ^a ROSAT All 
Sky Survey (RASS) mm a h ao (2000±400) BL Lac 
o6jeKaTa. 
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